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I. INTRODUCTION. 

The problem of electromagnetic currents in deuteron, 
especially of isoscalar nature, still cannot be considered 
as fully resolved despite of very numerous and longstand- 
ing efforts of many groups. As a good example we name 
here three topics from the field where the present ap- 
proaches all relying on the conventional nucleon-nucleon 
(NN) force models, failed to explain quantitatively the 
experimental deuteron data. These are: 

(i) the circular polarization P 1 of 7-quanta in radiative 
capture of spin-polarized neutrons in hydrogen at 
thermal energies |l|,|2|, viz- n+p — > rf + 7 reaction; 

(ii) the deuteron magnetic form factor B(Q 2 ) around 
the diffraction minimum at Q m in ~ 50 fm~ 2 0,0, 

HQ; 

(iii) the photon-induced polarization of the neutron 
from the d(^,n)p reaction at low energy Q. 

When going to few-nucleon electromagnetic observ- 
ables the disagreements of modern models for two-body 
currents with existing experimental data becomes even 
far more numerous |9( ■ The main part of the above dis- 
crepancies is related to the isoscalar magnetic currents 
which are now still strongly model dependent [S 0, H, 0, 
IToL ITU , and thus cannot be established uniquely. More- 
over, the existing theoretical approaches seems to employ 
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all the important types of currents, viz. onc-nucleon, 
two- and three-nucleon ones with many types of meson- 
exchange currents (p^j, unj, etc.) to remove the dis- 
crepancies. So, the various discrepancies which still re- 
mained could imply likely that some important contri- 
butions to the electromagnetic currents are still absent. 
These current components ignored until now are able to 
remove at least some of the disagreements observed to 
date in the deuteron- and few-nucleon electromagnetic 
observables. 

In the present paper we propose a model for new cur- 
rents of isoscalar nature, which still have not been dis- 
cussed up to date. It is demonstrated below that the new 
current removes the discrepancies enumerated above in 
(i) - (ii) that they make the theoretical framework for 
electromagnetic properties of the deuteron and the few- 
nucleon systems more consistent and thus more powerful. 

Hence, we shall discuss in Sect. [H] briefly the status of 
isoscalar magnetic currents in the deuteron and in the 
NN system at low e nerg ies (see more extended reviews 
in [H [H E E Q E3). The dibaryon model for the 
NN interaction is elaborated in Sect. Hill where we also 
include both EFT- and microscopic quark shell-model ap- 
proaches to substantiate our model. In Sect. lIVI the prop- 
erties of deuteron as emerging from the new force model 
are discussed in detail. We describe the short-range 
currents induced by the intermediate dibaryon compo- 
nent and the microscopic quark model formalism for the 
deuteron electromagnetic currents. Sect. E] is devoted 
to a general consideration of the isoscalar Ml- and E2- 
transition amplitudes and we present in detail the formal- 
ism for the isoscalar currents. In particular the formalism 
for the calculation of the np—^dry reaction, the form factor 
B(Q 2 ) and the magnetic moment p c i is discussed. The 



numerical results and their comparison to the respective 
experimental data are presented in Sect. IVII The sum- 
mary of the results obtained is given in Sect. EH In the 
Appendices some useful formulas and some details of the 
derivation have been collected. 



II. THE STATUS OF ISOSCALAR CURRENTS 
IN THE DEUTERON 

A consistent and correct description of isoscalar cur- 
rents in the deuteron and in the TVA^-system at very low 
energy seems to open a door toward quantitative studies 
of the short-range NN-iorce at low energies. Through 
this also toward to the non-perturbative QCD at low 
energies where there are no extra difficulties encoun- 
tered compared to investigations at high energies, such 
as a necessity for relativistic treatment, incorporation of 
many inelastic processes etc. The problem has been first 
posed by Breit and Rustgi 0] many years ago. Main 
part of the modern efforts to understand the problem 
quantitatively comes from effective field theory (EFT) 
treatment in high-orders (N 3 LO, N LO etc.) and also 
from the model treatment of short-range iVW-interaction 
within phcnomcnological iVA^-potentials. The point is 
that the main contribution to the amplitude of processes 
like n+p — > d+j comes from the isovector Ml current. 
This two-body current is dominated by a long-range tt- 
meson exchange. This phenomenon has been named in 
literature as a chiral filter hypothesis ^?|]- So, the isovec- 
tor Ml-transitions are "protected by the chiral filter" and 
does not manifest any sensitivity to short-range interac- 
tions. 

In the sharp contrast to this, two-body isoscalar cur- 
rent is not protected by the chiral filter and thus depends 
sensitively on the short-range interactions and short- 
range current contributions. This renders such currents 
to be highly model dependent. As a result, within the 
framework of EFT-approach the treatment of isoscalar 
currents requires some high orders of ChPT-expansion 
(N 3 LO and N 4 LO), which, in turn, introduces some extra 
parameters On the other hand, in a more 

conventional MEC-treatment (see e.g. 0, 0, 0|) the 
quantitative description of isoscalar currents includes, in 
addition to the usual p7r7-contribution, also the relativis- 
tic corrections dependent on the NN interaction model. 
All these contributions depend crucially upon the meson- 
nucleon form factor cutoffs A nNN , ^ttNA, A pNN , A^nn 
etc. and also upon the electromagnetic form factors of 
the intermediate mesons 10]. It is important to empha- 
size that the above cutoff values chosen for the two-body 
current models are generally not the same as in the in- 
put of NN and 3N potential models |2J]j and rather often 
are simply fitted to the electromagnetic observables mea- 
sured in experiments. The clear evidence to the strongly 
enhanced values for the cutoffs A^nn, A^jva, A^n etc. 
chosen in traditional iVTV-models like the Bonn poten- 
tial model can be seen in the results of relativistic cal- 



culations for the deuteron structure functions A(q 2 ) and 
B(q 2 ), and the deuteron magnetic moment 10]. It is 
very likely that a serious over estimation of the relativis- 
tic contributions found in Ref. ^(j indicates to the too 
high momentum cutoffs in the meson-nucleon vertices. 

Thus, one can summarize that the quantitative treat- 
ment of isoscalar currents in all existing approaches is 
related to the rather strong model dependence, in con- 
trast to the isovector Ml current. For example, the ex- 
perimental value for the circular polarization P 1 of 7- 
quanta emitted in the np — > process measured some 
times ago Q is P^ xp ~ (-1.5 ± 0.3) x 10~ 3 and is un- 
derestimated by all existing NN- and two-body current 
models |21l . However, despite of some inner inconsisten- 
cies in the two-body current models mentioned above it 
occurred to be feasible to describe many deuteron elec- 
tromagnetic observables like static characteristics, charge 
and quadrupole form factors, the np — > g?7 cross sections 
at low energies etc., rather accurately. But the descrip- 
tion of other electromagnetic observables, in particularly 
those given above in the items (i) - (iii), meets quite 
serious difficulties. 

In particular, there is a long standing puzzle tightly in- 
terrelated to the deuteron isoscalar current is a behavior 
of the deuteron magnetic form factor B(Q 2 ) in the area 
around the diffraction minimum, Q 2 ~ 45 - 55 fm~ 2 . 
There is a huge literature devoted to calculations of the 
deuteron form factors. However, the most recent fully rel- 
ativistic treatment of the -B(Q 2 )-behavior [fj has clearly 
demonstrated that the existing two-body current mod- 
els have missed some important contributions. Missing 
the same contributions very likely is the reason for the 
visible disagreement of vector and tensor analyzing pow- 
ers in p-d and n-d radiative capture at very low energies. 
Thus, it is evident, there is a number of mutually in- 
terconnected effects in the deuteron and the few-nucleon 
systems where one needs a new isoscalar current. How- 
ever, the conventional two-body current models still have 
no appropriate candidates for this. 

In this paper we propose an alternative mechanism 
for isoscalar current generation in the two-nucleon sys- 
tem based on the dibaryon model for NN interactions 
at intermediate and short ranges developed recently by 
our group H3, 0> HH- Generally speaking, em- 
ployment of the dibaryon degree of freedom to describe 
quantitatively the electromagnetic deuteron properties at 
low energies is not new. As an example, we shall re- 
fer to numerous studies having appeared in 80-ies where 
authors tried to incorporate various types of six-quark 
bags to the deuteron wave function to explain some puz- 
zles observed in electromagnetic observables (see, e.g. 
Refs. |H IHEl IM E3). These early attempts have 
revealed the important role which the quark degrees of 
freedom plays in the short-range NN interaction and in 
the electromagnetic structure of deuteron (see, e.g., [3(j 
and references therein). The very important role played 
by a "hidden color" and the quark pair currents in the 
deuteron electromagnetic form factors has been estab- 



2 



lished in these works, in particular. However, these early 
attempts have resulted neither in a deeper understand- 
ing of the deuteron electromagnetic observables nor in 
significant improvement of the description of experimen- 
tal data. 

Very recently, however, in tight interrelation to a re- 
naissance of the dibaryon physics (which is related par- 
tially to a recent boom in pentaquark phy sics), m any 
new studies, both theoretic al ]3ll 132 , l33L 134 . l35l l36l l37| 
and experimental [sE l4(l||. have appeared in liter- 
ature. An EFT-approach which included intermediate 
dibaryons |3^,|3^,|33 focused just on description of elec- 
tromagnetic properties of the deuteron at low energies 
and momentum transfer. However, the dibaryons in that 
approach have been introduced more on formal basis than 
as a real physical degree of freedom (or an intermediate 
broad resonance), to improve the convergence of chiral 
perturbation series and to simplify the EFT scheme in 
that case. In our dibar yon approach, in contrast to the 
formal scheme [jnl US HE EDI- the dibaryons are con- 
sidered as an important new degree of freedom in the 
description of NN- and 3JV interactions in nuclei at in- 
termediate and short ranges I22T . I23T.I25]] . 

The dibaryon concept of the nuclear force at intermedi- 
ate and short ranges turned out very fruitful in a quanti- 
tative or qualitative description of numerous phenomena 
in hadron and nuclear physics. In particular, this model 
made it possible to explain without any free parameters 
the Coulomb displacement energy for A — 3 nuclei ( 3 H 
and 3 He) and all ground-state characteristics of these 
nuclei |4l|; with exception of the 3 H and 3 He magnetic 
moments where we included only single-nucleon currents. 
However, according to general principle of quantum me- 
chanics, any new degree of freedom must lead to new 
currents of diagonal and transitional type. On the other 
hand, in our approach the weight of the dibaryon com- 
ponent in deuteron does not exceed 2.5 - 3.5%. There- 
fore, just the transitional current from the NN to the 
dibaryon component can be more important than the di- 
agonal dibaryon current. Nevertheless, we have consid- 
ered in the present study the diagonal dibaryon current 
as well. 



III. THE DRESSED BAG MODEL FOR THE 
MEDIUM- AND SHORT-RANGE NN 
INTERACTION 

A few key points have been taken as the physical jus- 
tification of this dibaryon model. 

(i) Failure of the c-meson exchange to describe the 
strong intermediate-range attraction in the NN 
channel [43, Hj, |4J|. Instead of a strong attrac- 
tion required by the iVW-phase shifts, the exchange 
of two correlated pions with the S*-wave broad a- 
meson resonance in tt-tt scattering included to the 
mechanism, when treated consistently and with a 
rigor, leads to a strong medium-range repulsion and 



a weak long-range attraction. Thus, the crucial 
problem arises, how to explain the basic strong 
attraction between two nucleons at intermediate 
range. 

(ii) The heavy-meson exchange with the mass m pM ~ 
800 MeV corresponds in general to the inter- 
nucleon distances Vnn ~ 0.3 fm, i.e. it occurs 
deeply inside the overlap region between two nu- 
cleons. And thus, these exchanges should be con- 
sistently treated only with invoking a six-quark dy- 
namics. The ignorance of this dynamics in conven- 
tional JViV-potential approaches (of the OBE type) 
leads to several serious inconsistencies in OBE- 
description of short-range NN force (see the dis- 
cussion in Refs. |22j, |23J, |45J ) and is also a reason for 
some problems with description of the short-range 
NN correlation in nuclei. 

(Hi) One more serious difficulty is related to the choice 
of different cutoff parameter values A^nn ; A^jva 
etc., when one set is used to describe the elastic 
JVJV scattering, another set to fit the pion produc- 
tion cross section in JVJV scattering and the me- 
son exchange contribution to the two-body cur- 
rent, and a third set to describe the three-body 
force 0, 152, . Using these different values for 
the identical vertices in the description of tightly 
interrelated quantities points toward an inadequate 
framework employed in the underlying model. 

The dibaryon model [231 |25| seemingly overcomes the 
above difficulties with a consistent description of the 
short-range NN force. Firstly, it employs a rather soft 
cutoff parameter A^nn (moreover, there is likely no se- 
rious sensitivity to these cutoffs in our model) . Secondly, 
the short-range NN interaction is described through 
a dressed dibaryon dynamics, so that the conventional 
heavy-meson exchange in i-channel plays no serious role 
in the short-range JVJV-dynamics. And thirdly, the short- 
range JViV-correlations in nuclei tested by high-energy 
probes can consistently be described as an interaction of 
the probes with the dibaryon as a whole, which may be 
associated with its inner excitations or de-excitations. 

In our approach the virtual dibaryon in the NN system 
are modeled through generation of a symmetric six-quark 
configuration s 6 [6] x dressed with a strengthened scalar 
(e.g. tr-meson) and other (e.g., n, p, ui etc.) fields. It can 
be schematically illustrated by Fig.^ The enhancement 




FIG. 1: Driving mechanism for the dressed dibaryon gener- 
ation used as an intermediate state in the short-range NN- 
interaction. 
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of the scalar field in the symmetric six-quark configura- 
tion s 6 [6] x implies some rearrangement of quark-gluon 
fields in the region, where two nucleons are totally over- 
lapping. The emergence of a strong scalar-isoscalar field 
in the six-quark bag induces automatically an isoscalar 
exchange current in the multi-quark system. 

The cr-meson loops originate mainly from "non- 
diagonal" transitions from the mixed-symmetry 2fuv- 
excited configurations |s 4 p 2 [/jsf][/cs]) to the unexcited 
fully symmetric configuration |s 6 [6x]) in the six-quark 
system with emission of a (virtual or real) cr-meson. 
In turn, the strongly attractive interaction between the 
scalar-isoscalar cr-field and the multi-quark bag results 
in an enhancement of the attractive very short-range di- 
quark correlations in the multi-quark system. Thus, as a 
net result of all these highly non-linear effects the mass of 
the intermediate dibaryon surrounded by the cr-field gets 
much lower as compared to the respective bare dibaryon 
(see the discussion below). 

Completely similar to the behavior of nucleonic system 
in a scalar field, where the scalar exchange may be viewed 
as a renormalization of the nucleon mass [46| 

m — ► to* = m + v s (0)p, 

where i> s (0) is a volume integral of the scalar exchange 
potential and p is the nucleon density, the constituent 
quark mass in multi-quark system is also reduced due to 
interaction with the scalar field: 

v s (0) glq/ml, 

where g aq is the cr-quark coupling constant; note v s (0) < 
. Thus, the dibaryon mass should be renormalized no- 
ticeably due to this scalar field mechanism. 

The effect of the strong attraction of the cr-field to 
the quark core and the resulting mass shift of the multi- 
quark system is illustrated easily by the anomalously low 
mass of the Roper resonance iV*(1440) with positive par- 
ity, which lies even lower than the lowest ne gat ive parity 
resonance iV*(1535). It is widely accepted [4|| that the 
Roper resonance has mainly the structure |iV + cr). How- 
ever, in the language of the simple quark shell model, 
the second positive parity level in the nucleon spectrum 
corresponds to 2Hlu excited three-quark configurations 
\sp 2 [21] x ) or \s 2 d - (s 2 2s)[21] x ). Thus, exactly the 
same mechanism as in our dibaryon model, i.e. the cr- 
meson emission from the 2fiw-excited state of the nu- 
cleon |sp 2 [21]x) — > |s 3 [3]jc + cr), leads to a generation 
of a strong cr-field and to a significant shift downward 
of the Roper resonance mass. Without this strong non- 
linear effect, the mass of the iV*(1440) would be ~ 500 
MeV higher than that found experimentally This very 
large shift gives some estimate for the magnitude of the 
attractive effects which appear in the interaction of the 
scalar field with fully symmetric multi-quark bag, or in 
other words, the effect of dressing by the cr-meson field. 
At very short NN distances the qq correlations become 
repulsive due to one-gluon exchange, and jointly with a 



strongly enhanced quark kinetic energy, this results in ef- 
fective short-range repulsion in the iViV-channel. In our 
approach this repulsive part of non-local NN interaction 
is modeled by a separable term An<£ag(r}^os( r ') with a 
large positive coupling constant Ao |22u23l | while the form 
factor (pos( r ) is the projection of the six-quark |s 6 [6]x) 
state onto the NN -channel: v?os(r) = (NN\s 6 [6] x )- 

Such combined mechanism lies fully beyond the per- 
turbative QCD, and we suggest it can be described 
phenomenologically by dressing the six-quark propaga- 
tor Ge q (E) with cr-meson loops 0, |25| . The resulting 
dressed-bag (DB) propagator Gdb(E) and the transition 
vortexes NN -> DB and DB -> NN treated within a 
microscopic 6q model |2^| lead automatically to a non- 
local (separable) energy-dependent short-range NN po- 
tential VMqN (r, r' ; E) 

V NqN = (NN\V Nq \DB( S 6 ))G DB (E)(DB{ S 6 )\V qN \NN) 
= <p(r)X(EMr'), (1) 

where the form factors ip(r) are deduced from the mi- 
croscopic 6c/ model while the coupling constant X(E) is 
determined by a loop integral with the cr-loop as shown 
schematically in Fig. ^ [23| . This loop integral can be 
conventionally parametrized in the following Pade form 



X(E) = A(0) 



E a + aE 
Eq-E- 



(2) 



Here the parameters A(0), Eo and a can be either calcu- 
lated from the microscopic 6q model or obtained from fits 
to the phase shifts of NN scattering in the low partial 
waves. This approach resulted in the Moscow- Tuebingen 
(MT) potential model of NN interaction p l l ^ . 

In a more general treatment recently developed in 
Refs. |3 HH the one-pole approximation for VNqN wa s 
obtained on the basis of a fully covariant EFT approach. 
In a more simple version of the model jHHa] the transi- 
tion operator V/vg which couples the nucleon-nucleon and 
dibaryon channels was calculated within a microscopic 
quark model with employment of the quark-cluster de- 
composition of the short-range N-N wave function. We 
consider here the NN — > DB transition in terms of this 
simple model, shown schematically in Fig. [21 where we 
assume that the coupling between the NN and DB chan- 
nels is realised on the quark level via a scalar exchange. 
This scalar interaction can be presented not only through 
the cr-meson exchange but also through a quark con- 
finement or another force including even the four-quark 
instanton-induced interaction of t'Hooft's type |47j . 

The transition operator VN q can be written in the form 



v Nq = ^2^ f g'iv a (n j ), 



(3) 



where u s (ry) is a scalar qq interaction. This opera- 
tor should be substituted into the transition matrix el- 
ement in Eq. 1(1)1. The particular form of the scalar op- 
erator and its origin are not significant here. The 
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specific mechanism of dressing can also be disregarded 
at this step since the dressing has already been taken 
into account in the propagator Gdb- The most impor- 
tant point is that the S-wave two-cluster state 3q + 3q, 
which represents the N-N system in the overlap region, 
should be very close in its symmetry to a superposi- 
tion of the excited six-quark configurations s A p 2 (see 



Refs. 0i C3 f° r detail) having non-trivial permuta- 
tion symmetries with the following Young tableaux in 
the coordinate- and color-spin spaces {/} = {[fx], [fas]} 
Q/ x ] = [6], [42], [/ cs ] = [42], [321], [2 s ], [31 s ], [21*]). For 
this reason, we can rewrite the effective NN interaction 
Eq. generated by the transition operator Eq. J2J in 
the following constrained form: 



V NqN (r,r';E) ~ ^{{iV(123)KJV(456)|}|«V{/}) (sV{f}\V Nq \s 6 [6] x ) 
ff 

x Gdb (E) (s 6 [6] x \V Nq \s 4 p 2 {f'}) (sV{f'}\ {|JV(123))|JV(456))}, (4) 
I 



where we leave the mixed-symmetry 2fiw-excited six- 
quark configuration s A p 2 only (but with all the possible 
Young tableaux {/}), instead of the total sum over all 
the excited six-quark configurations s 4 p 2 , s 2 p 4: , . . . , etc. 
Then, one can deduce from Eq. Q that the matrix ele- 
ment of the NN — > DB transition is proportional to the 
wavefunction of an excited nucleon-nucleon 2S*-state: 

{<jV(123)|<iV(456)|}| S y{/}) =C fi p 2S (r). (5) 

Here C/'s are purely algebraic coefficients and r = 
(ri + Y2 + r3 — r4 — rs — rg)/3 is the relative-motion 
coordinate of the two 3g-clusters. For simplicity, we 
use here the harmonic oscillator (h.o.) 2S'-state function 

V>2s(r) = 2r 3/2 Tr- 1 / 4 y /3/2(2r 2 /3r 2 -l)exp(-r 2 /2r 2 ) 
for the projection of the mixed-symmetry six-quark state 
onto the iViV-channel, which is characteristic of the quark 
shell model. The nucleon wave function N(ijk) in such an 
approach is described by the pure s 3 configuration of the 
CQM 

|JV(123)) = | S 3 [3] X {/ ST }> 

= ^(pi,p 2 )|[2 3 ] c [3] ST ), (6) 

where tpN(pi,p2) = Afexp [- ^i{p\ /2 + 2p§/3)] , the pa- 
rameter b is the scale parameter of the CQM, with b ~ 
0.5 - 0.6 fm; the relative coordinates arc p 1 = r\ — r 2 and 
P2 = ( r i + r 2)/2 - r 3 . 

Then, using the 2S function for the transition NN — > 
DB vertex in Eq. |JT)|. i.e. substituting ip(r) = <fi2s(?), 
we obtain in our simple ansatz Eq. <|3J) for the qq pair 
interaction u s (r.y) 

(NN(sV)\V Nq \DB(s 6 )) = g 2 (v) ^(r), 

and X(E) =gUv) 2 G DB (E) 1 (7) 

where (v) is a superposition (with the algebraic coeffi- 
cients Cf) of the quark shell- model transition matrix el- 
ements (s V | y^U SjU v * (y ) I s& ) ■ 

Note that we did not include the projector 
|s 6 [6] x ) (s 6 [6] x | onto the lowest six-quark configuration 



s 6 [6]x in the sum within Eq. Q, because this configu- 
ration, according to the most recent 6q calculations |50| 
lies well above 2.5 GeV and it is not touched by a strong 
renormalization, as does the mixed-symmetry state due 
to strong coupling to the er-meson field. Thus, the 
strong repulsive contribution from the bare intermedi- 
ate s 6 bag to the effective short-range NN interaction in 
the Moscow- Tuebingen (MT) model [2^, 123] is modeled 
by an orthogonality condition to the nodeless OS state: 

/ i>NN(r)vos{r)d 3 r = 0, (8) 

where (pog is a projection of the 6q-bag state onto the 
product of nucleon wavefunctions: 

Ws(r) =A/- f T 1 <7V(123)|(iV(456)|}| S 6 [6] x ). (9) 

Through the orthogonality constraint JSJ - © the sym- 
metric six-quark bag configuration is excluded from the 
NN Hilbert space, which prevents a possible double 
counting of the s 6 configuration. As a result, the total 
wavefunction of the two-nucleon system ^tot is defined in 
the direct sum of two Hilbert spaces Hnn © 'Hub- This 
direct sum can be conventionally represented by the two- 
line Fock column. For example, the deuteron state in the 
MT-model reads 

where the mixing angle 0N q is calculated on the basis of 
coupled channel equations with the transition operator 
VM q taken as a coupling potential [4lJ . Here the deuteron 
wavefunction in the ./VW-channel (or the iViV-component 
of the deuteron) takes the conventional form 

\d{NN)) = ^yl\ { {f) + ^-yfUf), (ii) 

r r 

where the S-wave component u(r) satisfies to the con- 
straint: 

/>oo 

/ u(r)cp os (r)dr = 0, (12) 
Jo 
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The deuteron iViV-component \d(NN)) and the dibaryon 
component \DB) are normalized individually to 1. This 
implies a standard normalization of the total wavefunc- 
tion 



(d\d) 



9 Nq (d(NN)\d{NN)) + sin 2 6 Nq (DB\DB) 



= 1 
(13) 



The dressed dibaryon propagator Gdb(E) can be rep- 
resented through the Dyson equation: 



G 



DB 



G 



(0) 
DB 



^DB^ GdB, 



(14) 



where Gp B is the bare dibaryon propagator and the £ is 
an eigenenergy which includes irreducible diagrams, i.e. 
those which do not include the intermediate free nucleon 
lines (but still can include NA or A A channels). Our 
DB-mode\ calculations kept only the leading er-loop 0j 
l2fij in the series for the E-kernel, because all other graphs 
(calculated within the six-quark shell model for the quark 
bag) corresponds to much higher masses (for the details 
see Refs. |24j,|25j). Moreover, the propagator for the bare 

six-quark bag Gp B corresponds to the pure six-quark bag 
states with Ofku and 2huj quanta for even partial waves 
and lhu for odd partial waves. 

When including the cr-meson loops only into the dress- 
ing mechanism, the dressed bag pr opagat or in the one- 
pole approximation takes the form |23 



[G 



DB\ LL > 



= / dk 



Bl(k)Bj,(k) 
E-E a (k) ■ 



(15) 



where B£(k) is the aDD form factor and E a (k) is the to- 
tal energy of the dressed bag (in the non-relativistic case 
E a (k) ~ m a + £ CT (k), with e cr (k) = k 2 /2m a + w CT (k) ~ 
m„„ + k 2 /2m r7 , and m aa — m rT m a /(m <7 + m a ) while 



= V? 



k 2 is the relativistic energy of cr-meson) . 



Thus the effective interaction in the TV TV-channel induced 
by the intermediate dressed dibaryon production takes 
the form in partial wave representation 



V NqN = £ Wi L ,(r,v';E) 



(16) 



JLL' 



with 



Wi L ,(r,r';E)^Y,^L M ^ XJ LL'(E)^ I (r'), 



A I 



X[ L ,(E)=j 2 [G DB }i LI , (17) 



where j 2 — gi(v) 2 (see Eqs. (@J - @ for comparison). 

It is easy to show [2^, EI that the weight of the DB 
in the total NN wavefunction in the LI/S 1 J- channel is 
proportional to the energy derivative 



All' j 



dXj L ,(E) 
dE ' 



(18) 



This derivation is in close analogy to a similar procedure 
for the weight of dressed particle state using the energy 
derivative of the respective polarization operator H(P 2 ) 
in the field-theory formalism. In particular, rece ntly we 
developed the fully covariant EFT-derivation |2J, [25J for 
the relativistic iViV-potential at intermediate and short 
ranges, similar to Eqs. l|15 jl -H18 jl : we also derived in the 
simplest one-state approximation a separable form for 
the relativistic iViV-interaction in channels 1 S'o and 3 Si- 
3 Di, which fits almost perfectly the respective NN phase 
shifts for the large energy interval - 1000 MeV. In con- 
trast to other potentials, e.g. the purely phenomenolog- 
ical Graz separable NN potential which includes a few 
dozens adjustable parameters for a lesser energy interval 
this high-quality fit has been performed using only four 
parameters in the singlet 1 5'o-channel and a few more for 
the coupled-channel case 3 Si — 3 D\ 

The separable form of the short-range NN interac- 
tion given in Eqs. I jlfijl - ljl7|) can be clearly compared 
to the contact terms in the effective field theory (EFT) 
approach (pionless) where all the high-energy physics is 
parametrized via some contact terms (see Fig. the 
left graph) which cannot be calculated within that low- 
energy approach and must be either parametrized some- 
how phenomenologically or fitted to the data p4l |25| . 
On the other hand, our short-range mechanism (shown 
schematically in Fig. EH right) gives just the energy de- 
pendence for such contact terms. 



IV. DEUTERON STRUCTURE AND 
DIBARYON- INDUCED SHORT-RANGE 
CURRENTS 

Similar to the general description of the iViV-system 
given in Sect. ITTTl the total deuteron wavefunction in the 
DB-model has the form of the Fock column Eq. I jlOjl 
with at least two components, the NN and the dibaryon 
dressed with tr-field — a dressed bag (DB) j23,|23j. The 
DB-component \DB > in the second line of Eq. (|10(l can 
be presented in the graphic form as a superposition of 
a bare dibaryon (six-quark configuration s 6 [6x]) and an 
infinite series of cr-meson loops coupled to the s 6 -quark 
core propagator. Thus, this component includes both the 
bare and dressed parts. 

Taking further the simple Pade-form J2J) of the energy 
dependent factor \ J LL , (E) in Eq. ijT?|) and by calculating 
the energy derivative fl%)) one gets easily the weight of 
the dressed dibaryon component in the deuteron, which 
turns out to be 1 ~ 0.025 - 0.035 for different versions 
of the model. It is very interesting that this weight of the 
dibaryon component derived from the energy dependence 
of the cr-loop diagram is rather close to the weight of non- 
nucleonic components (e.g. AA etc.) in the deuteron 
deduced within many phenomenological models. Table |I] 
gives the summary of the static characteristics of the 
deuteron found with this dibaryon model (2^, [^lj . 

The parameters of the model obtained from fit to the 
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TABLE I: Deuteron properties in different models. 





E d 


Pd 




Qd 




A s 


D/S 


Model 


(MeV) 


(%) 


(fm) 


(fm 2 ) 


{(in) 


(fm" 1 / 2 ) 




RSC 


2.22461 


6.47 


1.957 


0.280 


0.8429 


0.8776 


0.0262 


AV18 


2.2245 


5.76 


1.967 


0.270 


0.8521 


0.8850 


0.0256 


Bonn 2001 


2.22458 


4.85 


1.966 


0.270 


0.8521 


0.8846 


0.0256 


DB {NN only) 2.22454 


5.42 


2.004 


0.286 


0.8489 


0.9031 


0.0259 


DB {NN+Qq) 


2.22454 


5.22 


1.972 


0.275 


0.8548 


0.8864 


0.0264 


Experiment 


2.22454(9)° 




1.971(2) b 


0.2859(3) c 


0.857406(1; 


) d 0.8846(4) e 0.0264 / 




"Ref. [52], 




"Ref. [53], c Ref. [54] and 


[55], d Ref. [56], 


e Ref. [57J, 


•'Ref. [58]. 



phase shifts of elastic NN scattering in the 3 5i - 3 Di 
channel have the following values 

A(0) = -385.89 MeV, E a = 855.29 MeV, 



a = -0.025, r = 0.38113 fm, {b=d-r ). (19) 

We emphasize here again, that after fitting the phase 
shifts our approach does not have any free or adjustable 
parameters left for the description of the deuteron prop- 
erties. From Table H] one can see that the predictions 
for the basic deuteron observables found in the above 
dibaryon model is generally even in better agreement 
with respective experimental data than those for the 
modern iVTV-potentials calculated conventionally. 

The modeling of the NN — > DB transition by scalar 
exchanges between quarks makes it possible to consider 
the "contact" jNN — * DB vertices (Fig. gj) in terms 
of CQM with the minimal electromagnetic interaction of 
the constituent quarks, i.e. with the quark current 



(20) 



where q = p' — p, Si = \+\tz and F q {q 2 ) is a form fac- 
tor of the constituent quark which can only show itself 
at intermediate momentum transfer Q 2 > 1 GeV 2 /c 2 . It 
implies that the constituent quark is an extended object 
and has obtained its own electromagnetic form factor, 
e.g. a monopole one F q {Q 2 ) = 1/(1 + Q 2 /A 2 ), where the 
parameter A 9 is expected to be set by the chiral symme- 
try scale A x ~ 4nf n ~ 1 GeV. 

For definiteness, we consider the diagrams (a) and (c) 
depicted in Fig.0] In our model with the scalar exchanges 
and the CQM current l|20|) these contact terms are equiv- 
alent to the sum of Feynman diagrams depicted in Fig. 
These diagrams describe the two-particle currents in the 
six-quark system. The diagram (e) in Fig.JSJgives rise to 
an additional (i.e. the 7-induced) contribution to the 
transition s A p 2 — ► s 6 as compared to the mechanism 
shown in Fig. [21 Here we demonstrate that within the 
above model for the short-range NN interaction the min- 
imal quark-photon coupling j q {q)A fJi {q) leads to a non- 
additive two-nucleon current which does not vanish in 
the low-energy limit qo, |q| — > 0. In this limit, only the 



contribution of the diagram Fig. |Sfe) vanishes because of 
orthogonality of the s 6 configuration to the quark-cluster 
states in the N-N channel (i.e. to the configurations s A p 2 
and the other ones). By contrast, the total contribution 
of the diagrams {a)-{d) and (/) in Fig.[3]does not vanish 
at qo ~~ * 0. In each pair of diagrams, (a) and (6), de- 
picted in Fig. El the singular terms ~ l/<7o are mutually 
canceled, but the remainder, proportional to the momen- 
tum k between the ?-th and j-th quarks and also to the 
scalar qq- interaction potential i> s (k 2 ), does not vanish in 
the limit qo — ► because of non-vanishing the matrix 
element (s 4 p 2 \v s \s 6 ) (see below). 

Now we pass to the actual calculations of such dia- 
gram contributions. It should be stressed here that the 
current diagrams in Fig. corresponds just to the 'non- 
diagonal' (transition) electromagnetic current which cou- 
ples two different channels, i.e. the proper NN- and DB- 
channels. These channel wavefunctions enter the transi- 
tion matrix elements with a proper own normalization 
because any current associated to the DB-st&te is "nor- 
malized" to the weight of the Di?-component. Among 
other things this makes it possible to avoid any double 
counting (symmetry properties of quark configurations 
are discussed in Appendix^which argue strongly against 
the repeated contribution in detail) . When the spin part 
of the i-th quark current (Eq. I20fl is taken into account 
only and a low energetic Ml-photon is generated, one 
can write the following Feynman amplitudes M- for the 
diagrams depicted in Fig. |5fa) and (6) 



ieg 2 v s {k 2 ^ 



M 



2m q + fa— A z.„nv „ 



-e^rq,e^*\u{ Pi ), 



v (b) 2m 



-2pi • q 

u(p'j) <j ija^q^e^ 



2p>-q 



+ 2m :^/ g 3 -^4 A) *}»fa-). (21) 

where ej^ is a space-like photon polarization vector 
£ ai(a)* _ {0,6^-**} satisfying the transversality condition 
q-e (A) * = at A = ± 1." 

It is easy to verify that the singular terms 
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a^ 1 ' q v e]i '*m q /(p' i ■ q) and a^ v q v s\t '*m q /{—pi ■ q) cancel approximation qo/m q <C 1 the three-dimensional opera- 
each other in sum MHs + MH b -> in the limit qo 0. As tor 
a result, we obtain from Eq. (|21|l in the non-relativistic 



= tp w - - {«•«> - « - ^) + «•«> (- ^ x « - ^) } 

x (2^) 3 <5 3 (p l +p J -p' i -p;.-q) (22) 



defined on non-relativistic quark wavefunctions of the 
CQM. This operator describes the transition from the 
NN channel to the 6q-bag with emission of a Ml 7- 
quantum, i.e. a "contact" NN — * DB + 7 interaction, 
schematically shown in Fig.^Ja). 

Now we can calculate the effective contact vertex 
NN 4^ NNy (see Fig. BJi) on the basis of the quark 
operator V^q-f in terms of the quark-microscopic version 
of the DB model. Recall that in our model the diagrams 
in Fig. taken without electromagnetic insertions are 
simply the pairwise q-q interaction: 



which describes the transition from the NN- to the 
6g-bag channel (see Fig. |2J. This observation points 
toward the proper solution of the problem of contact 
NN <^> NNy (further on we use the notation "NqNy" 
for brevity) interaction in our approach. Namely, we cal- 
culate the non-local iVqA^-interaction operator in the 
NN Hubert space Vjv 9 Ar 7 (r, r') by the same way as 
the non-local iVgiV-interaction operator Vjv ? Ar(r, r') in 
Eqs. Q -0. 



3 6 
2 \ ^ \ ^ 



Vn q 9 s 2^2^[vs(K) + v s( k j)] We obtain finally ( see Appendix B for detail) the 

3 3 , , NqNy (contact) term searched for (as the sum of two 

x (2tt) S ( Pi + Pj - Pi - Pj - q), (23) gra phs, (a) and (6), in Fig. HJ| 



') = 



eZ 
2M N 



x q 



x q 



i 1 -/ ,^ dip 2 s(r) H E ') , 
x 1 -Ji(W 2 ) — j ¥>2S{r ) + (f2s{r)— ji{qr /2) 



■e^GUq 2 ) 
X(E)1 



dr 2M N 



'2M N q- 



d<f2s(r') 
dr' 



(24) 



where Gf f (0) = \i v + \x n and G\ t (0) = \i v — fi n (the ori- 
gin of the nucleon form factors Gff and G v in the quark 
model results of Eq. (|24|l type is discussed in Appendix 
C). Our basic expression for the transition dibaryon 
current still does not take into account possible effects 
which must affect the predictions of our model (viz. the 
relativistic effects and quark boost contributions which 
should be essential at Q 2 ~ 1 GeV 2 /c 2 @, and other 
contact terms with inclusion of pseudo-scalar and vector- 
meson exchanges 0,0 etc.). So, to account of all these 
effects we renormalize our contact NN NNy vertex 
using some renormalization factor Z in Eq. I|24l) . It is felt 
that the value Z»l± 0.3 is reasonable since a precision 



of 10 - 30% is typical for standard quark model evalua- 
tions of the hadron magnetic moments. We show below 
that when choosing a reasonable value for a single free 
constant Z — 0.7 the contact term given in Eq. (|2~4"|) leads 
to a considerable improvement in description of isoscalar 
magnetic properties of the deuteron. 

One can get a general expression for the electro- 
magnetic current in the NN system, and also in the 
deuteron, starting with the quark current 112011 which has 
already been used for finding the contact NqNy vertex 
PH>. When the total deut eron wavefunction in the Fock- 
column form Eq. I|10l) is considered, the diagonal matrix 
element of the quark current (|20|) can be represented as 
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(d\J23i(lH X) (Q)\d) = co S 2 8 Nq (d(NN)\J^ + V}f%\d(NN)) 

6 6 

+ sin 3 6 st (DB\^j?eW\DB) + 2co S 6 NqS in6 Nq (d(NN)\Y,J^i X) \DB). (25) 



The last two terms in Eq. I|25|) are nothing else but con- 
tributions of the graphs shown in Figs. |S£f) and (e) re- 
spectively. 

It is worth to summarize here our main findings. 
Within our two-component NN interaction model, the 
minimal substitution leads basically to two different two- 
particle currents (in addition to the single-nucleon cur- 
rent written in the first term of Eq. H25l) ~). viz. the 
transition (NN — ► DB) contact term as given by Eq. 1|24|) 
and the standard six-quark bag current given by the sec- 
ond and third terms in Eq. (|25p. 

In the nucleon sector we further replace the quark- 
model current of the nucleon J^qm with the standard 
representation of in terms of the phenomenological 
form factors given in Appendix C by Eqs. (|C1|) and l|C4() . 
However, the two-body current in the last two terms of 
Eq. (|25|l (which gives only a small correction to the single- 
nucleon current Jfi) is calculated here on the basis of the 
constituent quark model (see Appendix C for details). 

V. ISOSCALAR Ml AND E2 TRANSITION 
AMPLITUDES 

The effective electromagnetic operator of the isoscalar 
current Vj^ N derived in the previous section is defined, 
by construction, in a Hubert space of the NN component 
of the whole two-component system. Thus, it should be 
bracketed between the initial and final states just in the 
A^Af-channel. So, we look here at the application of this 
new current operator to the three observables: (i) radia- 
tive capture ri + p — > <i + 7 of spin-polarized neutrons by 
hydrogen; (ii) the deuteron magnetic form factor B(Q 2 ) 
in the region of its diffraction minimum; and (Hi) the very 
tiny correction to the magnetic moment of deuteron. 

In the radiative capture process, our main interest lies 
in the calculation of the circular polarization of 7-quanta 
emitted after capture at thermal energies. It includes 
both Ml and E2 isoscalar transitions. We can contrast 
for this process the "contact" isovector and isoscalar 
transitions, where the corresponding 7r-exchange term 
has a long range and corresponds to the isovector transi- 
tion while the scalar-exchange term relates to the short- 
range a-, 2-7T- or glueball-exchange between quarks in 



both nucleons (or to the instanton-induced interaction 
as well) and corresponds to the isoscalar transition. The 
above 7r-meson isovector current contributes to the total 
isovector amplitude for the So Si transition, which 
is generally large, and thus this term does not affect 
strongly the P 7 -value which is governed just by an inter- 
ference between isovector Ml and isoscalar M1+E2 am- 
plitudes. The main point here is that the single-nucleon 
isoscalar transition is strongly suppressed due to orthog- 
onality of the initial and final radial wavefunctions. In 
this case, the small isoscalar contribution to the P 1 and 
to an angular asymmetry of the photons can be of crucial 
importance due to their interference with the large isovec- 
tor amplitude 0, 0, . The isoscalar Ml current can 
be also very important for the deuteron magnetic form 
factor B(Q 2 ) in the area where the contribution of single- 
nucleon current almost vanishes. So, this new isoscalar 
current can affect essentially the behavior of B(Q 2 ) near 
its minimum. 

To fix uniquely the relative signs of partial transition 
amplitudes (and for a meaningful comparison between 
predictions of different models) we use in all our calcula- 
tions a common expansion of the photon plane wave into 
vector spherical harmonics (see, e.g. Refs. [6(3, |^) and 
a standard choice |62( for the phase factors of Clebsch- 
Gordon coefficients and spherical functions. This choice 
fixes the sign of the E2 amplitude uniquely. The prob- 
lem with the relative sign of E2 amplitude would arise 
if one calculates the E2 and Ml amplitudes separately 
(see, e.g., the detailed discussion of the i?2-sign problem 
in Refs. HI). 

This general formalism, common for two different elec- 
tromagnetic processes, has been used in the present work 
jointly with our new NN-lovce model to estimate a non- 
additive two-body current contribution. 



A. General consideration 

Let us start here with the single-nucleon current. The 
expansion of the circularly polarized 7-quanta emission 
(with A = ± 1) operator into electric and magnetic mul- 
tipoles takes the form [§3, 0| 
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J N ■ e^*(q) e-^' 2 = -= £ VMa + l)H) ! Jn • {j;(W2)^,A (0 
* z=i 



■iA 



I 



j i+1 (gr/2)>7 i+1 , A (f) 



(26) 



in which we employ for vector spherical harmonics 



Y 



E 



(i(A-x)lx|jA)Y I); ^(f) ?<*>($)). (27) 



with the spherical polarizations 



e (±1) (?) = T 



V2 ' 



K0)ivrt = 



(?) = 9, 



(28) 



which are determined in the reference frame related to 
the photon with the quantization axis Zq directed along 
the photon momentum q, zq = q. 

At thermal energies of neutrons, which corresponds to 
the velocity v = 2200 m/s, one can neglect the elec- 
tric dipole transition as the initial P-wave is strongly 
depressed. In Eq. (|26f) only two lowest multipoles, Ml 
and E2, remain, because the initial 3 D!-wave is admix- 
tured to the basic 3 iSi-wave by the strong tensor force; we 
emphasize here that our model includes also a significant 
short-range tensor force originated from the intermedi- 
ate dibaryon in addition to the conventional OPE tensor 
force. The operator of single-nucleon transverse current 



eJ N — 



1 



2m n 



£p±£li(-2iV T ) 



+e[(ji p a p + /i„<7„) x 2V r ]| (29) 

contains the V T , a transverse component of a gradi- 
ent which operates on the initial np wavefunction and 



a gradient V r = ^S/ r /2 operating only on the plane wave 
e -ig r/2^ associated with the momentum of emitted pho- 
ton q. Inserting ej£ into Eq. lj2"6l one gets, after some 
algebra, the following representation for the transition 
amplitude np — > in an arbitrary coordinate frame 
XYZ 



A=±l 



(30) 



Here the photon emission angles are given in the refer- 
ence frame XYZ where the axis Z is chosen conveniently 
along the polarization vector of incident neutron. Then 
the quantum numbers MM X' are projections onto the 
quantization axis Z of the initial and final spin of the np 
system and of the photon total angular momentum I re- 
spectively, while the photon helicity A = ±1 is defined, 
as above, in its own reference frame XqYoZo. In corre- 
spondence with this definition of the reference frame in 
Eqs. H26fl and l|30[l the matrix elements for the Ml- and 
E2-transitions are calculated for the fixed values A = ±1 
of the photon helicity, but in an arbitrary reference frame 
XYZ in which the initial-state wavefunctions of np scat- 
tering are given 



J 



- ^oir, p n )Y 00 (r) 



(iA p ±A n |00)xA p XA,A/,o 



r 



~ 3 Si(r,p„):Or) + -"D^p^lAr) (31) 
r r 



and the final-state wavefunction reads 

*lr(*0 = ^3&W + ^3&(f), 



(32) 



r r 
where y^(f) = £ (LM L SM s \JM)Y LML (r)XsM S . 

M L , JV/ S 

The states 1 S (r,p n ), 3 Si(r,p n ) and 3 Z? 1 (r,p„) are fixed 
by the dynamics of the np-interaction, i.e. the contin- 



uum 5*-waves are normalized by the respective scattering 
lengths 

1 S (r,p n ) -> r - a t , 3 Si(r,p n ) -> r - a s , 

Pn^O, (33) 

while the 3 D 1 (r, p n )-wave normalization is fixed by the 
3 Si — 3 Di tensor mixing in the initial state. 
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After elementary but lengthy calculations, one gets the 
following formulas for the matrix elements at the right 



hand side of Eq. (jSUJ 



Ml 



(A) 
MM' 



2m AT 



+(fj, n +flp) 



jo(qr/2)i 



(A)* . U P 



+l(jo(qr/2)+h(qr/2)) ■ (34) 



The first term in the curly brackets corresponds to the 
isovector Ml transition 1 S {NN) -> d( 3 5i) while the 
remaining two terms describe the isoscalar transitions 
in the coupled 3 Si - 3 Di channels generated by the 
spin-dependent and convection currents respectively. In 



contrast to this, the E2-transition amplitude is purely 
isoscalar, although it consists of two terms, convection 
(first) and spin-dependent (second) ones, similarly to the 
Ml transition, 



E2%,= [ dV*fr(f)^j- [\\j {qr/2)+n{qr/2)]eW*-- 2 L 



2 (»n+H P )j2(qr/2)V^Y 2 * 2>x (r) ■ -E-L-H. (35) 



The total np — > dj reaction cross section for unpolar- 
ized neutrons can be expressed through the respective 
amplitudes - l|55^l in the following way 



^ot m " 1=1 « 1 



\M1 



E E 

MM' X=±l 

(A) |2 
MM', 7=0 1 



4-7T 



\Ml { Xr =1 ? 



\E2 



(A) 
MM' 



(36) 



For the further calculations one can use the well-known 
properties of the Wigner 2?-functions, which gives actu- 
ally the angular behavior of the interference term be- 



tween Ml and E2 amplitudes. For example, in radiative 
capture of spin-polarized neutrons by spin-polarized pro- 
tons, n + p — » d + 7, the angular anisotropy for emission 
of 7-quanta in respect to the spin polarization axis of the 
initial nucleons can be described. 

Let us consider now the asymmetry in circular polar- 
ization of 7-quanta on the basis of Eqs. I|3()J1 and 134JI - 
H36J) . as measured in the experiment y|. The differential 
cross section for circularly-polarized 7-quanta emission in 
forward direction (i.e. along the spin polarization A„ of 
incident neutron) can be written in terms of the helicity 
amplitudes J31J - ESJ 



(Ta(A„) = —a\q\ 

Pn 4m 



2~o EE 

N Z X v M> 



M 



.|00)M1^ M , 



=1 + (iA p |A„|lM) 



Ml 



(A) 

MM', 7=0 



E2 



(A) 
A7A7' 



(37) 



For sake of brevity, the differential cross section respect to the neutron polarization vector in case when 
da(X ni 9 = 0)/dfl for 7-quanta emission at zero angle in the neutron spin projection onto the quantization axis 
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equals to X n is denoted as a\(\ n ). In Eq. (|3"?|) we have 
omitted the Wigner 2?-functions depicted in Eq. 1)30(1 be- 
cause the respective sums over A' are reduced, at 6 = 0, 
to a trivial factor 1. 

The differential cross section da(X n ,X p ;6)/dQ for the 
photon emission into an arbitrary angle 9 can be found 
using the simple transformation of Eq. ((37(1 by replacing 
the sum ^ Ea by Ea^i anc ^ by replacing the Wigner 
P-functions on the r.h.s. of Eq. 1(37(1 . 



B. The ft + p — > d + 7 reaction 

Using the general formulas (|36|l and 1(3711 for the helic- 
ity dependent cross sections one can find the circular po- 
larization P 7 (A„) and angular anisotropy r\ for the fixed 
initial values of X n (or A„,A P ) 



l D\ channels. Thus we get eventually 



Pj(X n ) — 



&unpol 



Acta (A n ) 

hJ2 E ^a(a„) 

A„ A=±l 



?y(A„, A p ) = 

da ( X„ , A p , 6 = § ) /dCl - da ( A„ , A p , 6 = 0) /dCl 
da(X n , X p , 6= § )/dft + da(X n , X p , 6 = 0)/dfl 



(38) 



(39) 



The Ml and E2 amplitudes which contribute to the 
cross sections 1)36(1 and (13711 are given in Appendix D 
in their exact form. It is important to stress, that the 
dependence of the Ml- and E2-transition matrix elements 
upon the momentum transfer q in Eqs. 1(D1(1 - ((D6(l is 
rather weak at low energies and becomes quite significant 
only for e-d scattering in the region of moderate and high 
momenta transfer (see below). This means that when 
applying formulas (|D1(I - (|D6(1 to the np — ► d-f cross 
section at the thermal energies the integrals I2 can be 
neglected while jo(qr/2) in the integrand of Iq can be 
replaced by unity. 

As a result, eventually the expression for P 7 can be 
presented via more simple "reduced" matrix elements 



Ml I= i = (fi p - fi n )I a (u, 1 S ), 
M1 I=0 = (ll p + ii n ) 



I (u, 3 S 1 )--I (w. 3 D 1 ) 



+ -I (w 1 3 D 1 ), 



E2 



1=0 



-W^i)- (40) 



with I (f,Z) = r o °° f(r)Z(r,p n )dr, where Z(r,p n ) can 
be any of the scattering wavefunctions in 1 S , 3 Si or 



P 7 (A„) = (-l) 1 /2+Anp. 



7' 



P 1 = 2Re 



Mlc 
Mil 



E2 Q 
Mil 



(41) 



in which the factor (— iy/ 2 + x ^ in front of P 7 reflects 
only that dependence on A„, which is deduced from the 
Clebsch-Gordon coefficient (iA p iA„|00) = (-l) 1 /^"^ 
in the first term of Eq. (|37|l . Moreover, since in the limit 
p n — > all ratios of the matrix elements in Eqs. (|Dljl and 
(|40|l become real, the symbol Re can be omitted here. 

Quite similar considerations of the angular anisotropy 
•q yield a formula quadratic with respect to the matrix 
elements in Eqs. 1(40(1 and bilinear on spin-polarizations 
of neutron and proton 0, 0] . 

In the literature there are calculations for P 1 with the 
RSC potential, but published results [2l| do not include 
any details and any patterns due to the interference of 
various Ml- and £ , 2-terms. Therefore, we make also a 
parallel calculation for the value of P 7 with the well- 
known Reid NN potential (version Reid 93 |13|)- Thus, 
the detailed comparison for all partial contributions be- 
tween our and the conventional RSC-potential model 
sheds light on the delicate balance of different isoscalar 
current components to the total value of P 7 . 



C. Deuteron magnetic form factor 

Usually the deuteron magnetic form factor includes a 
contribution of the transverse current fEa l29(l as a whole 
without explicit separation into Ml and E2 multipoles. 
However, in calculations of the deuteron magnetic form 
factor we still can employ the helicity amplitudes l(34|) 
and 1(35(1 derived above by summing and substituting the 
deuteron wavefunctions u and w instead of the wavefunc- 
tions 3 <Si and 3 D\ in the continuum. In this substitu- 
tion the isovector part (the first term of Eq. ((34(1 ) auto- 
matically vanishes while in the isoscalar part the electric 
charge e (in the convection current term) and the mag- 
netic moment (/x p + /i n ) in the spin-dependent term are 
replaced with their respective isoscalar counterparts, viz. 
the isoscalar electric and magnetic form factors of nu- 
cleon 



G s E (q 2 ) = G p E ( q 2 ) + G n E ( q 2 ), 
G s M (q 2 ) = Gl I (q 2 ) + G n M ( q 2 ), 



(42) 



As a result, the sum of Ml- and E2 contributions (see 
Eqs. I(D2() - ((D5(l in Appendix D) is transformed into 
the well known formula for the deuteron magnetic form 
factor 
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2 yZ-g 2 

3 2m n 



-G s E {q 2 ) j w 2 (r)[ ]a (qr/2)+ ]2 (qr/2)]dr 



u 2 (r) - \w 2 {r)^ j (qr/2) + ^ (u(r) + ^ 



h{qr/2) 



dr 



(43) 



where the factor y | accounts for the averaging of the 

amplitude squared over the spin projections. This gives 
the standard normalization of the deuteron magnetic 
form factor which leads to the conventional expression 
for the deuteron structure functions A and B 



eral formula l|34|) for the Ml-amplitude and modifies its 
spin-dependent part. The contribution of the contact 
dibaryon-induced interaction can be found if to replacing 
the isoscalar nucleon spin-current operator in Eq. I|34|l 



A(q 2 ) = [GUl 2 tf+ [GUl 2 )] 2 + [GUI 2 )]", 



B(q 2 ) = 2(1 + Vd ) [G d M {q 2 )\ 



Vd 



4mi 



(44) 



The cross section for elastic ed scattering is written as 

d(J e d OMott 



dn 



1 



m d 



A(q z ) + B{q z )tan z 



(45) 



The expression in the curly brackets of Eq. i|43|) evolves 
in the limit Q 2 — > and it goes to the well known formula 
for the deuteron magnetic moment 



Vd(NN) = n p + n n - -(/i p + fl n 

Pd -- 



D i 



w 2 {r)dr. (46) 



VI. THE COMPARISON WITH 
EXPERIMENTAL DATA 



leq 
2ttin 



A^n + Mp)jo(W2)e (A) * 



(47) 



with the respective spin-dependent operator for the 
dibaryon contact term 1)24(1 . Consequently, in the left 
bra-vector of the matrix element ((34(1 one needs to use 
the deuteron wavefunction instead of ty np . More- 
over, one can ignore in this case the energy dependence 
of the non-local potential V/v g .iv(r', r; E) in Eq. Q and 
substitute E — there instead of the real energy ex of 
thermal neutrons or the bound-state energy of deuteron 
Ed since the scale of the X(E) factor in Eqs. Q - (J5J) is 
of much larger and of the order Eq ~ 1 GeV. With these 
reasonable approximations one calculates first of all the 
isoscalar current contribution to the deuteron magnetic 
moment. 

A. The deuteron magnetic moment and the 
deuteron form factor 



When calculating the dibaryon- and quark contribu- 
tions to both physical processes one begins from gen- 



In our model the deuteron magnetic form factor takes 
the form: 



G d M (q 2 ) = y g 2M \cos 2 9 Nq G d M{NN) {q 2 ) + cos 2 B Nq ^ NqN F NqN (q 2 ) + sin 2 B Nq [i s eF s 6 (q 2 ) 

+2cos9 Nq sin9 Nq [j, s 6_ s 4 p 2F s 6_ s 4 p 2(q 2 )} , (48) 
I 

where the first term in the brackets represents the nu- NNj vertex ((24(1 
cleonic current contribution while the second one cor- 
responds to the isoscalar component of contact NN poo poo 

fi NqN F NqN (q 2 ) = G s M (q 2 )2Z / drdr'u{r)u(r') 



JO 



X(fi2s{r) 2M^q Jl{qr/2) ^^- (49) 
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Here, F NqN (0) = 1 by definition and thus the value fJ, NqN 
is equal to that of right-hand side integral at q = 0. 
The third and fourth terms in Ea. 1)48(1 represent the di- 
agonal and non-diagonal contributions of the 6q-core of 
the dibaryon, i.e. the bare dibaryon contribution. As 
is evident from Eq. I|C6)I of Appendix El the last term 
in Eq. f25|l vanishes at q = and thus does not con- 
tribute to the deuteron magnetic moment, while the ac- 
count of the second term in Eq. I|25[l results only in a 
minor renormalization of the deuteron magnetic moment 
(|46[l . As a result, the dressed bag gives a real contri- 
bution to the deuteron magnetic moment only due to 
contact NN NN^-vertex (I24|) . and this contribution 
is equal to 

A/j BB = cos 2 9 Nq n NqN . (50) 

In all the present calculations for the deuteron magnetic 
moment and the structure function B(Q 2 ) the published 
parameters of the Moscow- Tuebingen AA-model 

have been employed. These parameters allow to fit the 
NN phase shifts in the very large energy interval - 1000 
MeV. The mixing parameter 6 Nq can also be calculated 
with the MT model amounting to 

sin6 Nq = -0.13886. (51) 

The first term in Eq. l(4*5|) is calculated with Eq. (gBJ), 
while the third and fourth terms are calculated via 
Eqs. CH) - £U) of Appendix O The sum of these three 
terms to the deuteron magnetic moment amounts to 

/i d = cos 2 Nq fi d (NN) + sin 2 Nq fJ-d{6q) 

= 0.8489n.m., (52) 

as shown in Table 1. From the difference of this the- 
oretical prediction to the respective experimental value 
fi e d xp = 0.8574 n.m. one can find an admissible value for 
the contact-term contribution Eq. I(50|l , which in our case 
amounts A/i^ B ~ 0.01 n.m. The second term in Eq. 148(1 
is calculated with Eqs. 1491) - (|50|l and we employ the fixed 
values (O and (52) and Z = 1 in Eq. J22J). Thus, in this 
calculation for Afi BB we do not use any free parameters 
and reach a value 

A^ B = 0.0159, (53) 

which is in a very reasonable agreement with the above 
limitation. The resulting value for the deuteron magnetic 
moment ^ eor = 0.8648 n.m. overshoots the respective 
experimental value a little bit, but the remaining dis- 
agreement Afi = 0.0074 n.m. has decreased strongly. 

Now, in order to reproduce exactly the deuteron mag- 
netic moment, we fix the value of Z as follows: Z = 



0.7. We consider the accurate experimental value for 
the deuteron magnetic moment to give a stringent test 
for any new isoscalar current contribution. With the 
fixed renormalization constant Z = 0.7 we calculate the 
deuteron magnetic form factor and the circular polariza- 
tion P 7 . We follow this strategy in order to obtain a 
parameter-free estimate for the latter observables. 

Another restriction of the new isoscalar current is re- 
lated to the description of B(Q 2 ) near its minimum at 
Q 2 ~ 2 GeV 2 /c 2 corresponding to Q ~ 7 £m -1 . The 
position of the minimum depends crucially upon non- 
additive two-body contributions. In Fig. \7\ we display 
the results of our calculation for B(Q 2 ) based on Eq. l|4^|) 
and we compare them to the experimental data 0, HU- 
The dashed curve in Fig. represents the single- nucleon 
current contribution which is described by the term pro- 
portional to G d M ^ NN ^(q 2 ). The position of the minimum 
for this single-nucleon term appears noticeably shifted 
toward lower Q 2 values as compared to the experimental 
data 0. [3- Adding the conventional quark contribu- 
tion (dotted line in Fig. EJ) reduces this discrepancy due 
to the positive sign of the s 6 -bag contribution which is 
approximately compensated by the negative-sign inter- 
ference term between the nucleon and the bag contribu- 
tions. It is evident however from this consideration that 
one needs some positive contribution to reproduce the 
correct position of the minimum. 

In the model developed here, the contact term which is 
tightly related to the intermediate dibaryon production 
has just the necessary properties. Adding the contribu- 
tion of the DB contact term if^ in line with Eq. lf4*%|) 
results immediately in very good description for the 
deuteron magnetic form factor B(Q 2 ) as shown by the 
solid line in Fig. [7] Thus, a rather minor renormaliza- 
tion of the DB contact term by a factor of 0.7 makes 
it possible to describe quantitatively both the deuteron 
magnetic moment \i and the behavior of B(Q 2 ) in the 
large momentum transfer region Q 2 < 2,5 GeV 2 /c 2 . Fi- 
nally, by fixing this minor renormalization of the contact 
term the calculation the circular polarization P 7 will be 
parameter-free. 



B. The circular polarization of photons in reaction 

n + p — > d + 7 

The contribution of the dibaryon current to the 
isoscalar Ml transition 3 Si(NN) — > 3 Si(d) is calculated 
now in the same way. When the spin-dependent oper- 
ator 1471) in the matrix element (|34|l is replaced by the 
contact term <|24[) the Ml-transition amplitude for the 
circularly polarized 7-quanta emission is obtained as 
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AMI 



(A) _ 



MM' 



= Z / dV / d 3 rf d M ,(r') 



2M 



■'-IN 



x q 



■e W *G s M ( q 2 ) 



1 . f qr dtp 2S (r) X(E d ) , A(e r ) 1 ,qr' d<p 2S (r') \ , _ 



dr 2M 



N 



2M N q- 1 



dr 1 



When calculating P 7 this amplitude must be added to the 
single-nucleon current terms l|4U|l using the same renor- 
malization constant Z = 0.7. Similarly to the single- 
nucleon current, the integral (|54|l is calculated straight- 
forward by replacement of j\{qr /2) / 'q — > r/6 by q = E^. 
Also we substitute E = instead of E = er and E = Ed 
in the X(E) function in Eq. @. After this we get for 
the dibaryon induced current contribution an expres- 
sion analogous to Eq. (|D2|> with the respective "reduced" 
dibaryon matrix elements. The amplitude (|54|l found by 
this way together with the single-nucleon matrix elements 
(|40|l should be included to the final expression for P 7 



Tjtot 

1 



2Re 



AMI E2 



Mir 



Mli Mli MU 



(55) 



where the dibaryon induced current contribution is 

P 1 (DB)=2Re{^Y (56) 

The results of the numerical calculations within our 
model are presented in Table [H] together with a paral- 
lel calculation for P 7 with the conventional RSC NN- 
potential model in its modern version RSC93 [57j . Evi- 
dently the fully parameter free prediction of our dibaryon 
model for P 7 is in a first time in very good agreement with 
the respective experimental result. 



VII. SHORT DISCUSSION AND CONCLUSION 

In this paper we developed a model for the new electro- 
magnetic current in the deuteron and in the NN system 
in general. The new currents are based on the picture of 
short-range NN interaction via an intermediate dibaryon 
generation. The dibaryon represents a new degree of free- 
dom and according to a general principle of quantum 
theory this must inevitably lead to the respective new 
current (s). By applying the general recipe of minimal 
substitution to the Hamiltonian of the dibaryon model 
to derive the new current one gets automatically two dif- 
ferent contributions: diagonal and transitional ones. The 
diagonal current is associated mainly with the quark de- 
grees of freedom, and thus is proportional to the (small) 
weight of the dibaryon component in the deuteron. While 
the transitional current leads to a larger contribution to 
the deuteron electromagnetic properties, and likely also 
NN electromagnetic observables, especially of isoscalar 
nature. We studied three such electromagnetic charac- 
teristics: 



• the magnetic moment fid of the deuteron; 

• the magnetic form factor B(Q 2 ) in the region of its 
diffraction minimum; 

• the circular polarization P 7 of 7 quanta in radiative 
capture of spin-polarized neutrons by hydrogen. 

As for the prediction of the deuteron magnetic moment, 
the new isoscalar dibaryon current just fills perfectly the 
small gap which was found earlier (Afi ~ 0.010 n.m.) 
between prediction of the dibaryon TV TV-force model and 
experimental data (see Table QJ. With this tiny correc- 
tion the theoretical deuteron magnetic moment fid agrees 
excellently with its respective experimental value. 

In the present study we found that the minimal (gauge) 
substitution to the dibaryon Hamiltonian gives a strong 
positive contribution to the B{Q 2 ) behavior near the 
minimum region. Moreover, the parameter-free calcula- 
tion of the B(Q 2 ) in the new model gave already a very 
reasonable description for the deuteron magnetic form 
factor B(Q 2 ). A minor reduction of the dibaryon-7 ver- 
tex by a factor 0.7 results in an excellent agreement with 
the data both for fid and B(Q 2 ). 

After fixing all parameters of the new model, we calcu- 
lated the magnitude of the circular polarization of pho- 
tons inri + p—>d + "/ capture at thermal energy. This 
fully parameter-free calculation gave a result which is in 
a very close agreement with the existing experimental 
data [l| . It is important to remind the reader that many 
attempts were undertaken in the past; see e.g. the re- 
view of M. Rho 0| where one can find the references 
to earlier works and a good discussion of all difficulties 
encountered in theoretical predictions of P 7 . Thus, this 
longstanding P 7 -puzzle seems now to be solved. 

Here it is useful to discuss briefly the comparison be- 
tween the present model predictions and some other 
current models, both microscopic and phenomenologi- 
cal ones. Very detailed six-quark microscopic calcula- 
tions in Ref. [2!j have revealed that the quark exchange 
currents cannot give any quantitative agreement with 
deuteron data neither for the magnetic nor for the charge 
form factors, B(Q 2 ) and A(Q 2 ), respectively. More- 
over, when calculating the quark-exchange current cor- 
rections to the m agne tic and quadrupole deuteron mo- 
ments the authors [2£| have found some (although small) 
underestimation for fid but strong over estimation for Qd- 
These disagreements with the respective experimental re- 
sults have demonstrated that the incorporation of a bare 
six-quark contribution only cannot fill the gap between 
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TABLE II: Circular polarization of 7 quanta in the n + p — > d + 7 reaction 





P 7 (M1) 


P 7 (£2) 


P 7 (NN) 


P-y(DB) 


ptot 

7 


Model 


■10~ 3 


■10 -3 


•10 -3 


•10~ 3 


■10~ 3 


Reid 93 


-1.761 


0.699 


-1.062 


0. 


-1.062 


Moscow- Tuebingen 


-1.791 


0.657 


-1.134 


-0.261 


-1.395 


Experiment [lj 










-1.5±0.3 



the impulse approximation (plus the traditional MEC) 
results and the experiment, at least for the Ml- and 
752-isoscalar transitions. On the other hand, the dress- 
ing procedure for six-quark bag has been shown in the 
present work to lead inevitably to new short-range cur- 
rents. These dibaryon induced currents should replace 
the conventional two-body meson-exchange currents at 
short TViV-distances when two interacting nucleons are 
overlapping strongly in which case their meson clouds 
will fuse into one common cloud of a dibaryon. 

The new dibaryon currents proposed and studied in 
this paper must contribute also to many other elec- 
tromagnetic properties which could not be explained 
with the conventional AW-models before, e.g. the 7- 
induced polarization of nucleons at photo-disintegration 
of deuteron at low energies, d(7, n)p, and also the electro- 
disintegration of deuteron, d(e,e'p)n at high momentum 
transfer [g4| . The particular interest for the new isoscalar 
current rests in numerous studies of (e, e'pp) and ("f,pp) 
processes at intermediate energies. It is worth to re- 
member here that the theoretical interpretation of such 
processes, measured experimentally at various kinematic 
conditions, failed to explain the accurate experimental 
findings (see e.g. Ref. |Si|). Very likely these processes 
include some contribution of two-body isoscalar currents 
as well. 

Simultaneously, any success in such consistent inter- 
pretation of the data will support strongly the underlying 
dibaryon model for the short range NN interaction. 
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APPENDIX A: AVOIDING DOUBLE-COUNTING 

It is worth to add here some important comments 
about a large difference between "diagonal" and "non- 
diagonal" (transition) currents in the NN system. It 
is important that the similar consideration of the t- 
channel two-body currents in the TV TV system (with re- 
placement quarks by nucleons in Fig.[SJ) within the frame- 
work of a traditional NN model does not lead to new 
currents because the scalar exchange in the i-channel 
is already included somehow into the AW wavefunction 
^nn, and thus the one-nucleon current matrix element 
(^nn\Jn\^nn) also includes among the others the dia- 
grams shown in Fig. [B] 

In contrast to this (one-channel) problem, we are deal- 
ing here with a two- (or more) channel problem and our 
non-diagonal current operator VN ql given in Eq. de- 
scribes electromagnetic induced transition between two 
different channels, i.e. between the proper NN- and 
dressed bag components. In such transitions there hap- 
pens a strong rearrangement of the spin-isospin structure 
of the total six-quark wavefunction, and thus this transi- 
tion between two components is associated with the real 
current of magnetic type, which corresponds to the quark 
spin (isospin) flip. It can be illustrated clearly by consid- 
ering a transition |s 4 p 2 [42]x) — ► |s 6 [6]x) even with the 
same value of the total spin S (and isospin T) in the left 
(bra) and right (ket) vectors in the current matrix ele- 
ment: (s 6 [6]x,ST\V Nql \s 4 p 2 [A2} x ,ST). Recall also that 
the configuration s e in our model is fully absent in the 
initial and final NN states because of the orthogonality 
condition ljl2|l and the projector appearing in the NN 
channel. 

Let us assume that the spin and isospin in the initial 
(and final) N Af-channel and in the intermediate dibaryon 
state have the values S = 1 and T = 0, which are asso- 
ciated with the spin and isospin Young tableaux for the 
six-quark system [42] 5 and [3 2 ]t respectively. As a re- 
sult, we get the following transition matrix element 

( S 6 [6]x[42] s [3 2 ] T ([/^ T ]),[2 3 ] c | 

x V Nqi | S V[42] X [42] S [3 2 ] T ([/ ST ]), [2 3 ]c) (Al) 

It is important to stress here that [fg T ] ^ [Jst], i-e. the 
spin-isospin structure in bra- and ket states are different, 
because the Pauli principle severely restricts the form of 
the spin-isospin tableau [f' ST ] on the left hand side state 
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reducing it to a single allowable state [42] s • [3 2 ]t — > 
[/st] = P 2 ]s'T- In contrast to it, on the right hand 
side state any spin-isospin Young tableau [42] 5 • [3 2 ]t — > 
[/st] = [51]sT + [41 2 ]sT+j3i]sT + [321] S T+[2 2 l 2 ]sT is 
permissible (see e.g. Refs. pty). 

Therefore even the static magnetic moments for left 
and right functions must be different. It is worth to re- 
member here that just the spin-isospin Young tableau 
[f]sT determines fully the nucleon magnetic moment //jv 
at S = 1/2 and T — 1/2, and the correct value /zv both 
for proton and neutron can be obtained only for the sym- 
metric ST-state [3]st satisfying the Pauli exclusion prin- 
ciple. Hence, in the transitions IjAljl we have real cur- 
rents describing the spin-isospin flip in the NN — ► DB 



transition, and thus, the current operator Vat 97 makes a 
non-trivial (two-body) contribution to the total current 
in the two-channel system NN + DB. 



APPENDIX B: DERIVATION OF THE CONTACT 
iViV O iViV7 TERM 

In deriving of Eq. (|24|l we start from Eq. (@J and sub- 
stitute the vertex TV (77 I122H in one of two matrix elements 



of the 



s p 



the r.h.s. of Eq. (@J. It gives 



s 6 transition < s 4 p 2 {f}\V Nq \s 6 [6] x > in 



J 



V NqNl ~ J2 [{AWlsVU}) {sV{f}\V Nqi \s & ) G DB {s 6 \V Nq \ S V{f}) (sV{f'}\NN) 
ff 

MNN\ S Y{f}() S Y{f}\V Nq \ S 6 )G DB ( S 6 \V Nql \ S Y{f'}()sY{f}\NN)} (Bl) 



(non-significant details are omitted here). The calcu- 
lation of this operator is performed with the fractional 
parentage coefficient technique (see Refs. 0, EE E3 f° r 
details) by factorizing the fixed pair ij of quarks with 
numbers i = 3 and j — 6. The space coordinates of 
this pair depend on proton and neutron center-of-mass 
coordinates, r„ and r„ respectively, and on the relative 



motion Jacobi coordinates p pl , p p2 , 
fined in Eq. |JBJ and below 



Pni and Pn2, as de- 



where ir P 2 and n n2 are the relative momenta conjugated 
to p p2 and p n2 respectively, P = Y!Li Pi- 



r 3 = R + r p - 2p 2 /3, r 6 = R + r„ - 2p n2 /3, 



1 6 

R =e5>- ( B2 ) 



>=i 

In the cm. frame R = and r p = — r„ = r/2. The 
conjugated momenta P3 and P6 reads 

p 3 = P + 2(p p - p„)/3 - 7Tp 2 , 

p 6 =P-2(pp-p„)/3-ir„2, (B3) 



The 



calculation of the Fourier 
• • / 4^ ea; P(*Pi • r x + + ip 6 



(2^)3' 

elements 



transform 
• r 6 ) for 



of operators Vjv g7 and Vjvo on 



(27T) 

the matrix 

the r.h.s. of Eq. IjBlfl and substitution of the integral 
/ exp{— i(p3 + P6 — p 3 — Pg — q) • x}<ix for the <5-function 
(27r) 3 (5 3 (p3 + P6 — p' 3 — p 6 — q) lead, after the simple 
but tedious mathematics, to the following expression for 
the r-dependent part of the NN — > D£? + 7 vertex in 
Eq. ijBTJl 



2M, 



-F N (q 2 ) g 2 a (v) je-^'/^-l-Vr • [e (A) x er p ] - _L Vr ■ [e W X «x n ]j ^(r). (B4) 



The terms with factors ^— '■[o'i x q] on the r.h.s. of with the interpretation of the vertex matrix element ii 
Eq. (|22|l are canceled after the summation in the Eq. J7J. 
Eq. 1B1|) . Now the expression 1|B4|) should replace the 
function <p(r)ip2s(r) in the separable potential in line 
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In Eq. (|B4|> we use the standard formulas of CQM for 
the nucleon form factor 

F N (q 2 ) = F q (q 2 ) [ |7MPl,P2))|V^W3 d 3 pid 3 p2(B5) 



and for the magnetic momentum of nucleon 
(N(3q), [21] S [21] T : [3] ST | 



= Mw^jv (B6) 

7 2M N ' V ; 



N = n, p, with /i p = 3, /i n = -2. The gradients in Eq. 1B4|) 
originate from the proton and neutron momenta k p and 
k n which appear in the momentum representation of the 
vertex (|B4|I as a result of substitution of Eqs. I|B3|) 
for the quark momenta P3 and P6- In the coordinate 
representation, these momenta transform into gradients 
— iV rp = — «V r /2 and — iVr„ = *V r /2 respectively. The 
nucleon mass M N in Eq. I|B6(1 is a result of the substitu- 
tion of M N for the value 3m q which appears in the r.h.s. 
of Eq. (|B1|) resulting from CQM algebra. 
Making use of the expressions 



ar 

j a (qr/2) ± iV^h(qr/2) £ Y lm (q)Y* m {f) + ..., (B7) 



and adding the contributions of two graphs (a) and (b) 
depicted in Fig.0]we obtain finally Eq. (21} for the NqN^f 
(contact) term. 



APPENDIX C: NUCLEON AND DIBARYON 
ELECTROMAGNETIC MATRIX ELEMENTS 

The general formula for the nucleon electromagnetic 
current 



are not distinguished seriously from their experimentally 
measured counterparts used in Eq. (|C1() : 



H P = 2.79, it* = -1.91, 

F 2 (Q 2 ) « F dlp (Q 2 ) 



(C4) 



It is important to stress here that the quark-model mag- 
netic form factor, as can be seen in Fig. El resembles quite 
precisely the respective nucleon form factor found with 
the standard dipole fit, in the region Q 2 ~ 1 - 2 GeV 2 /c 2 , 
i.e. near the region of minimum for the deuteron mag- 
netic form factor B(Q 2 ) — see the respective curve in 

Fig. Elfor the choice of b = \f\ r o — 0.47 fin. 

Thus, in the nucleon sector, we replace the quark- 
model current J^qm with the standard representation 
of the nucleon current J% given by Eqs. ljCl|) and IIC4I) . 
and the non-additive two-body current (which gives only 
a small correction to the single-nucleon current J%) only 
is calculated on the basis of the constituent quark model 
(CQM). 

In particular, the CQM technique is used for calcula- 
tion of the last two terms in Eq. Ij25(l . These terms are 
contributions of the graphs shown in Figs. E{f) and (e) 
respectively. 



(i) For the graph in Fig.^f), the diagonal matrix ele- 
ment (DB\ Yl^Ui Ji £ ^ \D B) i s reduced here to the 
matrix element with the s 6 bag-like wavefunctions. 
This term is found by the same technique as for 
the nucleon matrix element calculation in Eqs. (|C2|) 
and (|C3|) . As a result, for the transverse current 
component (A = ± 1, — {0,e^}) one gets 



eJ%=e N rFf(q 2 ) + i 



2Ms 



<r^q»F 2 N (q 2 ), 



e« = e(i + ^) (CI) 



can easily be derived after averaging the quark current 
(|20|l over the nucleon wavefunction: 



^ QM = {N(123)\^(q)\N(123)) 



a .e/i N Q M e N 
6 " 7 +l 2M N ° q 



in which the quark-model predictions 



F QM (q 2 ), (C2) 



f^pQM — 3, 



2M 



(Q 2 ) = F q (Q 2 )e- Q2 " 2/6 (C3) 



(DB\ y £tf E M\DB) = 

6 

(s 6 ,S=l, T=0\J2tfeP>\8*,S = l, T = 0) 



(a p +a n ) ■ e 



-» s N F s e(Q 2 ), (C5) 



with FAQ 2 ) = F q {Q 2 )e-^l 2 \ fi% - Mp + Atn . 

(u) With the graph in Fig. 0V), the transition ma- 
trix element NN — > DB can be found similarly 
to Eq. jlf. i.e. by using the expansion over the 
six-quark shell-model states with restriction by the 
most important low-lying states 
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(d(NN)\J2^e^\DB) = £ / M (r)(7V(123)|(iV(456)| S V{/}>^( S V{/}lE^4 A) l s6 



(<t p +<t») ■ {ul2S{NN)) J2 C f F s e_ s * p s (Q 2 ). (C6) 



Here F s e_ s i p 2(Q 2 ) 



5^ e -5QV/24 and < u |25(iViV) >= / °° u(r)p2 5 (r)dr. 
I 



APPENDIX D: DEUTERON A/1- AND £2- TRANSITION AMPLITUDES 



The Ml- and E2-amplitudes which are used for the cross section calculation in Eqs. I|36[l and (|37[1 take the following 
form (we present here them separately for the spin (s)- and convection (c) current components): 



m1 mm> 1 ( S ) = $m',-\<> M:0 (ji n - Hp) A [^^j J o(w, ^o) 



/ — ieq \ 



(Dl) 



Ml 



(A) 7=0 



(*) - <W-a (M« + Mp) [AV2(1M1 -A)|1(M-A)] 



x { IoK 3 Si) - i/ (^, 3 ^i) + -^(^W^i) + hziw^Di) } . (1)2 ) 



Mll^(c) = <W-> [AV2(1M1 -A)|1(A/-A)] | {/„(«;, 3 £>i) + I 2 («,, ^O} 



(D3) 



(1M2-A)|1(M-A) 



E2™.(<!) = ««.,„-, [AV2(1M1 -A)|1(M-A)] f^") | {/„(™, 3 -D,) + Hw.'D^} 

where Jo arL d I2 arc the following overlap integrals 

7 (/, Z) = / °° f(r)Z(r,p n )j (qr/2)dr, I 2 (f, Z) = J™ f(r)Z(r,p n )j 2 (qr/2)dr, 
hiuw^Sx^Dx) = i / °° [u(r) 3 D 1 (r,p n )+w(r) 3 S 1 (r,p n )] h(qr/2)dr. 



(D5) 



(D6) 



Here Z(r,p n ) can be any of the scattering wavefunctions in 1 S'o, 3 Si or 3 Di channels. 

For the sake of convenience, the factors which are equal to unity on modulo are separated out in square brackets: 



AV2(1M1 -A)|1(A/-A) = 1, 



(1M2-A)|1(M-A) = (-1) M . 



y/3 



When summing over the M and A induces, these terms play a role similar to the Kroneker delta. 
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FIG. 2: Graphical illustration for AW — + DB transition in 
terms of quark microscopic model. The double-dashed line 
denotes some scalar exchange which can induce the AW — > 
DB transition jHIH]. 




gS(r-r') <p(r)A,(E)<p(r') 

FIG. 3: Different ways for parameterizing the short-range NN 
interaction: in low-energy (pionless) EFT through the four- 
nucleon contact term (left) and in the dibaryon model (right) 
via s-channel intermediate dibaryon. 




(a) (b) (c) 



FIG. 4: Schematic representation of the new electromagnetic currents induced by intermediate dibaryon generation. 




(a) (b) (c) (d) 




(e) if) 



FIG. 5: The diagram series illustrating the contact VN ql terms contributing to the "non-diagonal" current. 
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p — R — T^R' p i — fi ~ p i p i ^ p i 

1 Ak i S + S 

Pj — - -,y ij — ■ -p- \) • |-i>- p , > ■ -iy 

>lq >lq 

(») (b) 

FIG. 6: The part of diagrams shown in Fig. contributing to the non-diagonal dibaryon current with explicit notations of the 
kinematic variables. The double-dashed line denotes the scalar exchange. 




FIG. 7: The structure function B(Q) of elastic e-d scatter- 
ing. The impulse approximation (IA) result for the Moscow- 
Tuebingen potential model [2^| is shown by the dashed line. 
The sum of I A and the diagonal (s 6 — » s 6 ) and non-diagonal 
(s 6 — > s 4 p 2 ) bare 6q-contributions is shown by the dotted line. 
The total contribution of the IA+ bare dibaryon + DB (con- 
tact term) contribution is represented by the solid line. The 
data are taken from Q (Saclay), (SLAC) and |j| (JLab). 
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FIG. 8: Comparison of the simple quark model description for 
the nucleon magnetic form factor (with different scale values 
b) with those parametrized through the standard dipolc fit. 
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